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The vertex reconstruction toolkit RAVE has been extended by an option for the inclu- 
sion of kinematic constraints, and embedded into the ILD analysis framework Marlin. 
The new tools have been tested with an exemplary reconstruction of WW and ZZ de- 
cays. The presented results show the improvements achieved in precision of the fitted 
masses, and demonstrate the usage and functionality of the toolkit. 

1 Software elements 

1.1 The Rave library 

The Rave library [5] was created with the aim of avoiding repeated re-implementation of 
similar algorithms in every new reconstruction software, when specific modules can be for- 
mulated in an experiment-independent manner. This is certainly true for the modules doing 
the reconstruction of interaction vertices, e.g. of particle decays. 

The core of the library is an implementation of algorithms for geometric vertex recon- 
struction which were developed for the CMS experiment at the Large Hadron Collider, 
augmented by an additional simple and stable interface. The same approach was taken for 
this work: the algorithms doing vertex reconstruction with kinematic constraints [3] were 
extracted from the CMSSW framework and implemented in the library. 

1.2 The MarUnRave plug-in 

A plug- in has been developed for the ILD analysis framework Marlin [J , enabling users to 
access the geometric and kinematic vertex reconstruction capabilities of RAVE. Its name 
is MarlinRave, and it provides two new processors to the framework - RaveVertexing and 
RaveKinematics - one of which or both to be selected by the user. 

1.2.1 The RaveVertexing processor 

Geometric vertex reconstruction (by linear or adaptive filter algorithms) is accessed from 
Marlin through the RaveVertexing processor. The configuration of this processor requires an 
input collection containing objects of type EVENT: :Track. It takes a text string configuring 
the vertexing algorithm to be used together with its parameters, and the names of two 
output collections: one containing the fitted vertices, and the other containing the re-fitted 
("smoothed") tracks at their vertex. 

1.2.2 The RaveKinematics processor 

Its algorithm is based on a linear filter with Lagrangian multipliers for the constraints. 
Because the configuration interface of Marlin does not allow for arbitrary nesting, the 
RaveKinematics processor is designed to contain a flexible number of pre-defined kinematic 
decay topologies, thus facilitating the inclusion of new topologies as needed by the user. 
From the Marlin configuration file only the topology is chosen, whilst a text string allows 
the passing of specific parameters. 
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2 Reconstruction of W^W decays 



Generated Monte-Carlo masses 
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rave's kinematic capabilities were tested with a data sample consisting of 1376 events 
e+e" — > W^W~ 4 jets at y/s ~ 500 GeV, generated by Pythia [5] and processed for the 
LDC 2007 layout [6]; the "true" W masses are plotted in Figure 1(a) Events with a mass 
\mw — 80.32| > 6.5 GeV are eventually discarded, retaining 1008 events. 

To model the expected performance of 
jet reconstruction in ILC experiments, the 
following Gaussian errors have been applied 
to all four generated jets: 

as/E = 30%/Ve (1) 
= (70 sin(0) — 10 mrad (2) 

The direction resolution affects the 6 (polar 
angle) and 4> (azimuth) measurements. 

A known general problem with the use of 
jets as virtual measurements is due to their 
compositeness: only direction and energy, 
but not an absolute momentum, are well- 
defined. However, an initialization of the 
jet's momentum with its energy, followed by 
an appropriate inflation of the associated er- 
ror, gives satisfactory results. 

The constraints applied are those of en- 
ergy and 3-momentum conservation. They 
are explicitly written as (the subscripts 
identify the four jets): 
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(a) All W masses as generated by Pythia. 



All combinations 
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El + E2 + E3, 



Pi +P2+ P3 



Ei=V 
-P4 = 



(3) 
(4) 



So far no distinction was made between 
the four jets: the applied kinematic con- 
straints acted equally on all final states, 
and did not take into account which par- 
ticles they could have originated from. Now 
one has the task of associating the four jets 
into two pairs (there are 3 possible combi- 
nations), with each pair originating from a 
W boson decay. 

A trivial strategy is to assume that each 
combination is valid, and to calculate the 
two W masses accordingly. Figure 1(b) 
shows the results: each entry represents the 
W pair masses corresponding to one combi- 
nation of one event 
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(b) The W masses resulting from all possible 
combinations of jet association. 

Figure 1: Scatter plots of the W pair masses. 
The colours represent the number of entries in 
each bin. 



About 5% of the combinations yield unphysical mass values and do not enter the plot. 
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A comparison of Figure 1(b) with 1(a) suggests that, as a first step, dropping combina- 
tions where both W masses exceed a pre-defined cut hmit would significantly improve the 
performance of the association. This cut is chosen at a value of 130 GeV. 

For a second step, several options are 
feasible. A simple "equal-mass hypothesis" 
introduces in the kinematic fit, as an ad- 
ditional constraint, the requirement of the 
two fitted masses to be equal (because they 
both belong to W bosons). However, such 
a requirement would strongly favour com- 
binations along the 45° diagonal over those 
parallel to the axes, in contradiction to the 
true distribution of the W pair masses as 



shown in Figure l(a 



The most obvious improvement of such 
a hypothesis would be to model a selection 
requirement from two uncorrelated Breit- 
Wigner (BW) distributions@ But that is 
not possible in real- world scenarios, because 
the position parameter of the expected dis- 
tribution is exactly the (unknown) value to 
be determined by this kinematic vertex fit. 
Therefore, a compromise is suggested be- 
tween this idea and the "equal-mass hypoth- 
esis" above [8]: 
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Figure 2: Scatter plot of the W pair masses fit- 
ted with the similar- mass constraint; selection 
of the best jet association by the fit's pseudo- 
probability, applied to the reduced sample 
after the 130 GeV cut. 



The two W masses are known to not being equal, however, they are still picked from the 
same distribution. The likelihood of such a configuration is given by 



L{dc\dm, m) = L{dc\d,n) ■ L{mi{dm), m2{dm)\m) 



(5) 



Here, the first term holds the results of the kinematic vertex fit (parameters ac and 
covariances Vc). The second term represents the new "similar-mass hypothesis"; it models 
the distribution of the two masses dependent on the position parameter m, and it can be 
split into two symmetric terms if the masses are approximately uncorrelated. Although any 
p.d.f. can be used within the second term, a Gaussian model is certainly the easiest choice. 
Then, the objective function for the new fit can be written as 
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(6) 



where crt sa 9 GeV is a pre-defined scale (obtained by integration over the BW). 

The objective function M is minimized w.r.t. dm and fh and yields a pseudo-x^, the 
probability of which is used as the new second step selection criterion for finding the best 
jet association (out of the 3 possible combinations). 



More precisely: from o'(s) = J ds^ f dsBp{sA)p{sB) with p{si) oc IiW{Mw > ^i) El- 
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Entries 


1958 


Mean 


79.34 


RMS 


12.07 


Underflow 





Overflow 


2 


/ ndf 


52.25 / 42 


Constant 


133.7 ±4.0 


Mean 


79.58 ± 0.27 


Sigma 


11.4±0.2 



This strategy, after application 
of the first step (130 GeV cut 
mentioned above), results in 979 
events; the W pair masses are 
shown in Figure [2l 

The association performance of 
this method is characterized by 
type 1 and 2 errors of 1.6% and 
0.9%, respectively. 

Figure [3] shows the recon- 
structed W masses after kine- 
matic vertex fitting, including the 
similar-mass constraint |f| 

A Gaussian fit over the his- 
togram range 30 130 GeV yields Figure 3: Reconstructed W masses after kinematic ver- 
for mean and standard deviation: ^^^"^S with the similar-mass constraint. 
H{mw) = 79.58 ± 0.27 GeV and a{mw) ^ 11.4 ± 0.2 GeV. 

Scaled to 664k events expected in 4 years at the ILC (500 fb~^, 4-jet efficiency w 40%) 
[8], an accuracy of 0.014 GeV in determining the W mass may be achieved. 




3 Reconstruction of Z°Z° decays 

The other exemplary reconstruction was performed on a data sample of 994 Pythia events 
e+e" Z°Z° 4 jets at = 250 GeV, and processed for the LDCprime 2008 layout 
[9]. In contrast to the WW data of chapter 2, simulation included initial state radiation; 
at present, this cannot be accounted for in our kinematic reconstruction. Therefore, events 
with a total ZZ energy < 249 GeV are discarded, retaining 531 events. 

The ZZ data had the momenta of the jets correctly calculated by PFA, so they were used 
in this study. Apart from that, exactly the same strategics and methods as in chapter 2 have 
been applied. A comparison (not shown) of the reconstructed Z masses with and without 
kinematic constraints, respectively, reflects significant improvement. 
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Note that the four kinematic constraints (equ. 3 and 4) and the similar-mass constraint (equ. 6) 
together do not suffice to account for the redundant degrees of freedom. 
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